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ABSTRACT 

Fundamental mysteries remain regarding the physics of Type la supernovae (SNIa) 
and their stellar progenitors. We argue here that important clues to these questions may 
emerge by the identification of those SNIa that occur in extragalactic globular clusters- 
stellar systems with well defined ages and metallicities. We estimate an all-sky rate 
of ~ 0.1ry(L>/100 Mpc)^ yr~^ for SNIa in globular clusters within a distance -D, where 
r] is the rate enhancement per unit mass as a result of dynamical production channels 
that are inaccessible in the galactic field. If ry f» 2 — 10, as suggested by observations 
and theory, the combined efforts of accurate supernova astrometry and deep follow-up 
imaging should identify the >1% of nearby {D < 100 Mpc) SNIa that occur in globular 
clusters. 

Subject headings: galaxies: general — globular clusters: general — supernovae: general 



1. Introduction 

Many questions remain about the most fundamental aspects of Type la supernovae (SNIa 



including the triggering and hydrodynamics of t he explosions (e.g ., iHillebrandt &: Niemeveiil200d ) 



and the nature of their stell ar progenitors (e.g. , lYungelsonl |2005| ) . The increasing SNIa div ersity. 



with some very bright (e.g., iHowell et al.l 120061 ). some v ery faint (e.g. 



Kasliwal et al.l 120081 ). and 



some that do not follow the Phillips (1993) relation (e.g.. |Jha et al.ll2006l ) has energized the discus- 
sion of many possible formation scenarios. 



While at most a few percent of white dwarfs explode as SNIa (IPritchet et al.M2008l'l. there 
are only loose c onstraints on the specific binary evolution pathways (e.g., Ilben &: Tutukovlll984l : 
YungelsonI |2005| ) . There is a consensus that SNe la originate from thermonuclear ignition and 
burning of a C/0 white dwarf in a binary system. Yet it remains uncertain if the event is trig- 
gere d by accretion from a hydrogen-rich companion or from a merger with another white dwarf 
(see [Branch et al.lll995l ). referred to as the single-degenerate and double-degenerate scenarios, re- 
spectively. In the single-degenerate scenario, the main issues are the nature of the companion, and 
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if mass transfer can be sustained onto the white dwarf at the favorable rate in a sufficient number 



of systems to match the observed SNIa rate (jMaoa l2008l ) . In the double-degenerate scenario the 
primary issue is whether off-center carbon ignition can be avoided, as th is would transform the 
C/0 white dwarf i nto a Q/Ne/Mg w hite dwarf (e.g.. iNomoto &: Ibenlll985l ) rather than generate a 
SNIa (see however iRegos et al.ll2003l ). 



In fact evidence is amassing that multiple progenitor scenarios lead to SNIa. In particular, there 
is a clear disparity in the SNIa rates in late-type g alaxies with active sta r formation and elliptical 
galax i es that contain mos tly old (5-10 Gyr) stars ([Mannucci et al.l l2005l : IScannapieco k. Bildsten 



20051 : ISullivan et al.l 12009 1. This suggests that 10^ - 10^° yr can elapse betw een the birth of the 



progenitors and the explosi ons, a challenge for any single progenitor scenario (jDella Valle &: Livio 



1994 



Mannucci et al.ll2006l ). Multiple routes to SNIa may well lead to diversity in the explosive 



outcomes. For example, low-luminosity SNIa are most prevalent in early- type (i.e., E/SO) ga, laxies 



while the most luminous events occur only in star- forming galaxies (e.g.. iHamuy et al.lll996l ). 



We propose that the study of SNIa in globular clusters (hereafter, GCIa) may provide unique 
clues to understanding SNIa. Although globular clusters (GCs) and elliptical galaxies are both 
composed mainly of old stars, there are crucial differences between them. In a given GC, we are 
certain that all stars were born within 1 Gyr of eac h other, whereas e lliptical galaxies often show 
evidence for a substantial spread of stellar ages (e.g., iTrager et al.ll2000l ). Secondly, in an individual 
GC, the stellar metallicities are narrowly distributed, and thus the integrated metallicity is a good 
measure of the metallicity of any of t he constituent stars. Until recently, there wa s only one known 



exception {oj Cen) in the Milky Way (iFreeman &: Rodgerall975l : iBedin et al.ll2004l ) , but more recent 



work has found that a few massive (> 10^ Mq) GCs have helium-rich sub populations (e.g.. iPiotto 
20081 ). In any instance, GCs have metal licities low enough tha t a single GCIa detection would place 



ogn 

m 



stron g constraint on theoretical models (jKobavashi et al.lll998l : iHachisu et al 
20081 ). 



199S 



Piro &: Bildsten 



GCs are differentiated from the galactic field by their high stellar densities (often >10^ Mq pc~^) 
that trigger frequent close encounters between stars and binaries. Such encounters are responsible 
for the high incidence of exotic objects in GCs, includ ing X-ray binaries, rapidly spinning radio pul 



sars, blue stragglers and c ataclysmic variables (e.g., iHut et al. 



1992 



Sills et al 



2000 



19991 : 



Rasio et al 



Poole V &: Hutl |2006| ) . Dynamics will almos t certainly play an i i nportant role in the produc 



tion of GCIa progenitors (jShara &: Hurley 



2002 



Ivanova et al 



2006 



Rosswog et al.ll2008l ). likely 



increasing the GCIa rate per unit mass. 

We start in §2 by estimating the GCIa rate and discussing the possible dynamical enhance- 
ments. The observational challenges to finding a GCIa are discussed in §3, where we motivate 
that the maximum distance for such a search is 100 Mpc. We also explain the need for accurate 
astrometry of nearby SNe that will enable meaningful followup observations. We close in §4 by 
describing the implications of detecting even a single GCIa. 
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2. The Supernova Rate in Globular Clusters 



Globular clusters (GCs) have ~10^-10^ old (> 8 Gyr) stars inside a few parsecs, with a wide 
range of < O .SZc?) metallicities. All galaxies contain GCs, with total numbers scaling as ~100 GCs 
per 10^° (jAshman k Zepjll998l ). The common measure of the GC number density is the l^-band 
specific frequency, 



Sn = iVcclO' 



(1) 



where Ngc is the number of GCs, and Ady is the absolute l^-band magnitude of the galaxy 



([Harris van den Bergh 



ellipticals (jHarris 



1991 



1981 



Typical values of Sj\f are ~1 for spiral galaxies and ~2-5 for 
). Though standard, Sn is not the best choice for our purposes. We are 
most interested in the fraction of stellar mass in GCs, Fqq = MQc/Mg, where Mqc is the total 
mass of the GC system and Mg is the total stellar mass of the galaxy. Given the galactic stellar 
mass-to- light ratio Ty, Fqq is related to Sn by 



Fgc 



1.2 X lO^^SArmsT"^ 



(2) 



where is the mean GC mass in units of lO^M©, and we use Mv,q = 4.8 for the absolute 
magnitude of the Sun. Photometric studies of GC systems find ~ 2, and old elliptical galaxies 
have Ty = 3, so that Fqq ~ 2 x 10~^ for most ellipticals, while it is somewhat less for spirals. For 
some central dominant ellipticals at the centers of galaxy clusters, Sn can reach 10 (Harris et al. 
2009), corresponding to Fqq 10^^. 



Scannapieco &: BildstenI (j2005l ) and iMannucci et al.l (j2005l ) proposed that the SN la rate in 



a galaxy is the sum of two components, one proportional to the total stellar mass Mg, the other 
proportional to the star formation rate Mg. Such a model suggests that SNIa result from at least 
two evolutionary channels. For a particular galaxy, the two-component rate can be written as 



SNR(t) = AMg{t) + BMg{t) , 



(3) 



where A and B are constants. Using a well characterized sample of SNIa and host galaxies, 
Sullivan et al.l ( 20061 ) find basic agreement with eq. ([3]) and determine A = (5.3ibl.l) x 10~^^ 



yr 



and S = (3.9 ± 0.7) X lO""^ yr"^^© yr" 



The value of A is derived from SNIas in E/SO galaxies with no discernible star formation (i.e., 
-6 = 0). If these galaxies are truly old, with negligible star formation in the past ~5-10Gyr, then 
a reasonable first guess is that the same rate per unit mass also applies to GCs (the enhancement 
due to stellar dynamics is discussed below). Adopting the galactic value of A for GCs, the GCIa 
rate in a galaxy is AMgFQQ. An estimate of the local cosmic rate density of GCIas is obtained as 
follows. 

At low redshift, th e total -fT-band lumin osity density is jx — (5 ± 0.5) x 10^ Loi^- Mpc~^ 
{Hq = 70kms~^ Mpc~^: iKochanek et al.ll200ll ) and the fraction in E/SO galaxies is 40-50%. Given 
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the mean iiT-band stellar mass-to-light ratio Tx ~ 1, the contribution to the la rate density from 
old stellar populations is Ajk^ k- Here we let Tx = 1 ± 0.5, averaged over all galaxy morpho- 
logical types, where the uncertainty larg ely reflects a range of model assumptions, rather than 
measurement error (e.g., ICole et al.ll200ll ). We find A mTx ^ (2.7 ± 1.5) x 10~^ yr~^ Mpc~^. This 
is consistent with the recent measurement at z ~ 0.1 (jPildav et al.l 120081 ) . showing the dominance 
of the old stellar population for the local SNIa rate. The corresponding GCIa rate density is then 
^JK^ k{Fgc) ^ where (-Fgc) is the mean GC mass fraction over a large number of galaxies. If we 
adopt a plausible value of {Fqc) = 10~^ (see eq. [2]), we estimate a 

GC rate from mass alone ~ 3 x 10^^ yr^"*^ Mpc~^. (4) 

Under our given assumptions, we expect that the net uncertainty in this rate is a factor of ~2. 

However, it has been known for over 30 years that X-ray biii aries are >100 times more abundant 
per unit mass in GCs than in the disk (IClarklll975l : lKatall975l ). Dynamical interactions involving 
single stars and binaries occu r frequently in GCs and naturally account for this overabundance 
(e.g., iBildsten &: Delovd l2004l ). There are also excellent observational and theoretical arguments 
that dynamics shapes the GC populations of blue stragglers, milliseco nd pulsars, and cataclysmic 
variables fe.g.. lSills et"aDll999l : Easio et al.l[200ol : IPoolev k Hutlbood ). 



Recently, IShara &: Hurleyl (120021 ) and llvanova et al.l (l2006l ) explored the idea that SNIa pro- 
genitors can be formed by dynamical means in dense star clusters, leading to a mean enhancement 
of the GCIa rate per unit mass, r/. Shara & Hurley (2002) suggest that the number of WD- 
nondegenerate star binaries is similar in dense clusters relative to the field, although they found 
strong differences in the masses of the companion stars, which could be important in determining 
which binaries support stable accretion. The models described in Ivanova et al. (2006) suggest an 
enhancement of r/ ~ 1-7 for single-degenerate progenitors, where the range in ry reflects variation 
with metallicity and other parameters. 

In the double-degenerate case, Shara & Hurley (2002) showed that supra-Chandrasekhar WD- 
WD merger rate is over an orde r of magnitude higher in dense clusters than in the field. On the 
other hand, llvanova et al.l (|2006l ) suggest a more modest enhancement of r/ ~ 2. Based on recent 



studies of the prevalence of post classical novae supersoft sources in M31 GCs, iHenze et all (j2008l ) 
conclude that the nova rate in GCs may be as much as ten times higher than in an old field stellar 
population, and they suggest that GCIa may be detectable in future surveys. Overall, it seems 
conceivable that rj ^ 1-10 and that observed GCIa may help differentiate progenitors. 



3. Observational Considerations within 100 Mpc 

The maximum distance of interest is set by the need to find the underlying GC. GCs have 
a distribution of absolute magnitudes given by dN/dMy oc exp [-(My - Mvflf/2al], where the 
dispersion is ay ~ 1-1.5 in relatively bright galaxies with My g^i < — 20. Over a wide range of host 
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galax y properties, the mean is Myfl = —7.4 to within a few percent (e.g.. lHarrislll99ll : lJordan et al 



20071 ) and has been detected at the 100 Mpc distance of the Coma clu ster, where the tu rnover 



apparent magnitude (~27.6) is accessible by the Hubble Space Telescope (iHarris et al.ll2009l ). Ob- 
servations from the ground are presently limited to magnitudes of mv ^ 26, corresponding to the 
GC luminosity function turnover at distances of <50Mpc. 

For a given limiting apparent magnitude Vmax, the fraction of stellar mass in GCs brighter 
than Fmax is 



1 /•(Vmax-Mv)/V2cr^ 

fiV < Fmax) = ^ / dx e-"' , (5) 



where My = Myfi + DM - 0.4(ln 10)c7^, and DM is the distance modulus. When DM = 35 and 
Vmax = 26, as cry increases from 1 to 1.5, / increases from 0.25 to 0.62, while the fraction of globulars 
with V < Vtnax varies from 0.05 to 0.14. Even though the fraction of visible clusters may be small, 
the fraction of stellar mass contained within these clusters can be substantial, and typically exceeds 
50% when DM < 35 and ay takes its usual values of ~1.3-1.5. Hence, observations of a galaxy at 
100 Mpc for which the GC census is complete for clusters brighter than V = 26 finds ~15% of the 
globulars by number but ~ 60% of the stellar mass. This is an important point, since a dynamical 
enhancement in the GCIa rate may favor more massive clusters, making it of considerable interest 



to pursue GCIa within DM = 35 (jPoolev &: Hutll2006l ). At distances greater than f«100Mpc, it 



becomes extremely difficult to identify a significant number of GCs, and those detected in the 
outskirts of the galaxy will represent only a small fraction of the mass of the GC system. 

From eq. ^ and our discussion of the dynamical enhancement of the GCIa rate per unit mass, 
we estimate a local GCIa rate of «0.1?7 (D/lOO Mpc)^ yr~^. A plausible value of ~ 10 results 
in ~1 GCIa per year within 100 Mpc, which is about 1% of the total SNIa rate within 100 Mpc. 
What are the prospects for carrying out such a search? The typical SNIa within 100 Mpc would 
have a peak visual magnitude of mv ~ 16 — 17 and even the subluminous, 1991bg-like SNIa would 
be found at these distances in the upcoming wide angle (one-tenth of the sky) nearby SNe surveys 
(e.g., Palomar Transient Factory, SkyMapper and Pan-Starrsl). The SNIa yields from these new 
surveys, as well as the increasing numbers from targeted galaxy and cluster searches by LOSS- 
KAIT, CHASE, and ROTSE, makes the time right for explicit GC identification efforts. In some 
cases, prior HST or ground-based studies will have GC catalogs to cross-list locations. However, 
the typical case will require waiting until the SNIa has faded below the GC light. 

The few well-studied late-time SNIa light curves reach My ~ —7 after ~ 600 days (e.g., 
Sollerman et al. 2004; Lair et al. 2006), at which point their fade rates are 1.4 magnitudes per 100 
days in BVR (Sollerman et al. 2004) with colors oiV — —1 and B — 1/ ~ 0, much bluer than a 
GC. The / band decays more slowly (~ 1 mag in 100 days), again pointing to BVR for discovery. 
One possible way to first discern the presence of an underlying GC would be the detection of a 
modified BVR color evolution as the redder BVR colors of the GC (V — Rk, 0.4-1) begin to shine 
through. 
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After identifying a candidate GCIa, high angular resolution observations are required: (1) to 
confirm that the SNIa and GC lie along the same sightline, and (2) to minimize the likelihood 
that the SNIa occurred in the host galactic field, in front or behind the GC. If the GC and SNIa 
are found to overlap within a resolution element of diameter 6 (in arcsec) the probability that the 
SNIa occurred in the field is roughly Lq/Lqq = e, where Lqc is the GG luminosity and Lg is 
the luminosity in field stars within the resolution element. A definitive GGIa detection requires a 
small value of e, and a correspondingly small value of Lg. Since surface brightness, and thus Lg, 
generally falls with increasing galactocentric radius, the strongest GGIa candidates will be located 
far from the centers of their host galaxies. To illustrate this point more quantitatively, we assume a 
de V a.ucouleurs profile with a su rface brightness He = 19.5 mags/ arcsec^ at the half-light radius. Re 



(e.g.. iDiorgovski &: DavisI 119871 ). We further assume that the target GC has magnitude M = —7.5 
and that the GG is unresolved (the typical GG half-light diameter is ^C'l beyond lOMpc). With 
these assumptions, we find that the radius at which Lg/Lcc = e is given by 



Re 



34 3 , /D 

\ log — 

25 10 ^ V e 



(6) 



where Z?io = D/lOMpc. For e = 0.1 and Diq = ^/W (~31.6Mpc) we find R/Re ~3.4 when 
9 = 0.1 arcsec, and R/Re —0.78 when 6 = 0.02 arcsec. At a distance of 100 Mpc, the same two 
6 values give R/Re — 7.6 and ~ 2.3. However, at 1.0 arcsec resolution, R/Re — 14.7 even when 
DiQ = vTo. Since typical half-light radii are Re = 1-4 kpc, it is clear that ^O'.'l resolution is 
required to achieve modest R for small e, which can only be accomplished from space or with 
ground-based adaptive optics. Even then, the best GCIa candidates will be at i? > 10 kpc, which 
requires accurate astrometry of both the active SNe and the possible underlying GG. These limits 
highlight the value of extremely accurate astrometry for making GGIa measurements in the future 
with large ground-based telescopes. 



4. Implications of a Discovery 



A major open question is how a 10 Gyr old stellar population produces an appreciable SNIa 
rate, as this requires double-degenerate mergers or stable accretion from relatively low-mass donors, 
neither of which are currently-favored for SNIa production. While some elliptical galaxies show 
definit e signatures of relat ively recent low-level star formation within an otherwise very old system 



(e.g., iTrager et al.l l2000l ) . we can be confident that no new stars are forming in old GGs. A 



single, definitive GCIa detection would demonstrate that, in fact, SNIa do occur in truly old stellar 
systems. 

Secondly, any systematic trends in GCIa properties with metallicity contain information about 
the physics of the explosions. GCs are, with few exceptions, extremely uniform in their che mical 
compositions, and the same cannot be said of elliptical galaxies (e.g., iMehlert et al.l l2003l ) . Of 
course, the metallicity can vary a great deal between GCs, but the metallicity of an individual 
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cluster is readily determined from its photometric colors. The detection of a single [Fe/H] < 
— 1 GCIa would place stron g constraints on mode l s of stable white-dw arf accretion from a non- 
degenerate companion (e.g., iKobavashi et al.lll998l : iHachisu et al.lll999l ). and analysis of a handful 
of GCIa would strongly constrain models of SNIa lightcurves and pre-explosion simmering (e.g., 
Timmes et aPbood : IPiro Bildstenlboosl l. 



A potential complication is that the dense stellar environment of the GC may open up exotic 
paths to SNIa (e.g., Rosswog et al. 2008). We would hope that such an outc ome would be revealed 



in comparison betw een the GCIas and the field SNe, both in rates (e.g., IShara &: HurlevI 12002 : 
Ivanova et al.ll2006l ) and systematic properties, such as lightcurve shapes. If the rates are sufficiently 
high, GCIa may prove decisive in implicating nonstandard la progenitors such as double-degenerate 
mergers. 

Because it takes approximately two years for a typical la to fade to the luminosity of an average 
globular cluster, GCIa will require late-time observations. This is not common practice, although 
it would require minimal investment of telescope time for the closest supernovae, as the total rate 
is only ~ 10yr~^ within 30 Mpc. Accurate astrometry of these events is also critical to the later 
GC search. For the moment, we must appeal to the la archives. In the Sternberg catalog, there 
are 112 SNIa identified from 2005 to 2007 within z < 0.025 (7500 km/s). Of these, 34 are hosted by 
E/SO galaxies, and 6 in this subset are separated by more than 1' from the centers of their hosts. 
Only careful follow-up observations will tell if the first GCIa has already been detected. 
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